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ABSTRACT 

Recent development of galaxy surveys enables us to investigate the deep universe of high redshift. We quan- 
titatively present the physical information extractable from the observable correlation function in deep redshift 
space in a framework of the linear theory. The correlation function depends on the underlying power spectrum, 
velocity distortions, and the Alcock-Paczyriski (AP) effect. The underlying power spectrum is sensitive to the 
constituents of matters in the universe, the velocity distortions are sensitive to the galaxy bias as well as the 
amount of total matter, and the Alcock-Paczyhski effect is sensitive to the dark energy components. Measuring 
the dark energy by means of the baryonic feature in the correlation function is one of the most interesting ap- 
plications. We show that the "baryon ridge" in the correlation function serves as a statistically circular object 
in the AP effect. In order to sufficiently constrain the dark energy components, the redshift range of the galaxy 
survey should be as broad as possible. The survey area on the sky should be smaller at deep redshifts than at 
shallow redshifts to keep the number density as dense as possible. We illustrate an optimal survey design that 
are useful in cosmology. Assuming future redshift surveys of z ^ 3 which are within reach of the present-day 
technology, achievable error bounds on cosmological parameters are estimated by calculating the Fisher matrix. 
According to an illustrated design, the equation of state of dark energy can be constrained within ±5% error 
assuming that the bias is unknown and marginalized over. Even when all the other cosmological parameters 
should be simultaneously determined, the error bound for the equation of state is up to ±10%. 
Subject headings: cosmology: theory — galaxies: distances and redshifts — galaxy clustering — large-scale 
structure of universe — methods: statistical 



1. INTRODUCTION 

Observing galaxy clustering in deep redshift surveys is one 
of the most direct way to probe the evolution and struc- 
ture of our universe itself. Recent advances in observational 
technologies of deep surveys are spectacular, and enormous 
amount of information on deep universe must be flooded 
in near future. To understand the cosmological meaning of 
galaxy clustering in deep universe, theoretical analysis of the 
clustering of observed galaxies are necessary. 

The observable objects in the deep redshift surveys are 
fainter than the shallow surveys, and the sampling number 
density per comoving volume is smaller if we adopt selection 
criteria by apparent magnitude of galaxies. When one is in- 
terested in linear clustering regime, where complex nonlinear 
dynamics is not relevant, the mean separation of the objects 
should be at least of order 10/!~'Mpc. Recent advances of 
redshift surveys, such as the Deep Extragalactic Evolutionary 
Probe survey (DEEP; Davis et al. 2003), and the Sloan Dig- 
ital Sky Survey (SDSS: IYork e t al. 2000), etc. can map such 
fainter objects with sufficient qualities. By an efficient color- 
selection technique, the sampling number density can be im- 
proved. 

The deep redshift surveys, if the sampling density is large 
enough, have many advantages. First, the observable volume 
at deep redshifts is much larger than that at shallow redshifts, 
simply because the available volume is large in distant space. 
Second, the nonlinear clustering properties are less impor- 
tant at a fixed scale in deep surveys, because the universe is 
younger in deep universe and there is less time for the density 
fields to evolve nonlinearities. Third, while the observations 
of the nearby universe can map only the snapshot of the re- 
cent universe, the observations of the deep universe reveal the 
evolutionary properties of the universe. This is an indispens- 
able information to establish a consistent picture of the whole 



universe. 

The measurements of the cosmic microwave background, 
such as t he Wilkinson M icrowave Anisotropy Probe (WMAP) 
satellite fBennett et al.'"2003), provide the information of the 
universe at very high redshift, z ~ 1 100. On the other hand, 
the traditional galaxy redshift s urveys, such a s the Las Cam- 
panas Redshift Survey (LCRS. IShectman et al. 1996). IRAS 
Point Source Catalog Redshift Survev (PSCz. l'Saunders et all 
l2000h . and the AAT two-degree field galaxy redshift sur- 
vey (2dFGRS,|polless et al. 2001) reveal the recent universe, 
z ^ 0.2. Systematic investigations of the universe between 
these redshifts are undoubtedly promising ways to unveil the 
evolutionary feature of the universe. Especially, the galaxy 
redshift surveys of redshift range of z = 0.2-3.0 can be attain- 
able by present-day technology. The Kilo-Aperture Optical 
Spectrograph (KAOS) project' is proposed to survey galax- 
ies at those redshifts with sufficient density, using one of the 
existing Gemini telescopes. 

One of the most significant feature of the universe one can 
probe by such deep redshift surveys is the nature of the dark 
energy component of the universe. The dark energy compo- 
nent does not have important contributions at very high red- 
shifts. The effects of the dark energy can be probed by study- 
ing how the recent universe evolves with redshifts. One of 
the striking evidence of the dark energy was provided by the 
luminosity distance-redshift relation of type la supernovae 
jRiess et a^ ^, 1998: Perlmutteret_al, 199 9) at redshifts z ^ 1. 
The galaxy redshift surveys had not been considered to pro- 
vide constraints on the dark energy, since the galaxy surveys 
had been restricted to the low-redshift universe at z ^ 0.2. 
The ongoing and future galaxy surveys will break this limi- 
tation and the evolutionary features of the universe, includ- 
ing dark energy properties, will be searched by galaxy sur- 

' http://www.noao.edu/kaos/ 
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veys. A machinery of probing the dark energy by galaxy sur- 
yeys is based on the extende d Alcock-Paczyiiski (AP) effects 
jAlcock & Paczvnski' 'Ball inger. Peacock & HeavensI 

[1996; Matsubara & Suto 1996) on galaxy clustering in 
redshift space. Recently, several authors investigate 
the feasibility of detecting the nature of the dark en- 
ergy fro m the extended AP effec ts of the correlation 
function fMatsubara & Szalav" '2003") and the power spec- 
trum ( Blake & Glazebrook 2003; Seo & Eisenstein 200^ 
iHu & H aiman 2003; Under 200^. 

Theoretical predictions of the extended AP effect are easier 
when the power spectrum is used, because the dependences of 
the observed power spectrum with AP effect on cosmological 
parameters are not complicated once the dependence on the 
mass power spectrum is known. Including simultaneously the 
peculiar velocity effect ( Kaiser 1987) and the AP effect on 
the pow er spectrum requires only simple pro cedures to im- 
plement (i Ballinger. Peacock & HeavensI 1996 h. However, ob- 
servational determination of the reliable power spectrum re- 
quires a homogeneous sample, since the observed spectrum 
is a convolution of the survey geometry and the real density 
spectrum. On the other hand, the correlation function is more 
straightforward to observationally determine. Once the selec- 
tion function of the survey is known, the correlation function 
is obtained by simple pair countings even when the survey 
selection is inhomogeneous and the geometrical shape of the 
survey is complex. Although the theoretical calculation of the 
two-point correlation function in redshift space is not as ob- 
vious as that of the power spectrum, we now have an analyti- 
cal formula which make a numerical implementation very fast 
(|Matsubara 2000; Matsubara, Szalav, & Pope 2004). The last 
formula also contains wide-angle effects of the survey, which 
can not be naturally implemented in the power spectrum anal- 
ysis. 

The purpose of this paper is to provide a thorough analysis 
of the two-point correlation function in deep redshift space 
and its dependence on cosmological parameters. The cos- 
mological information contained in the correlation function 
is quantitatively investigated. As a result, we provide error 
forecasts of the cosmological parameters probed by the cor- 
relation function in a given redshift survey. Complementarity 
of samples at different redshifts is quantitatively addressed. 
The analysis in this paper provides useful information for the 
design of the future surveys at deep redshifts. 

This paper is organized as follows. In §2, a compact ver- 
sion of the analytic expression of the correlation function in 
arbitrarily deep redshift space, derived from a more general 
expression of the previous work is presented. Physical con- 
tents in the analytic formula are addressed in §3. In §4, series 
of deep redshift surveys which are within reach of the present- 
day technology are considered and expected error bounds are 
estimated by means of the Fisher information matrix. Optimal 
designs for the future redshift surveys are indicated. Conclu- 
sions and discussion are given in §5. Some basic notations are 
introduced in AppendixlAl Useful equations in calculating the 
Fisher matrix with marginalizations are given in AppendixlEl 

2. THE ANALYTIC EXPRESSION OF THE CORRELATION 
FUNCTION IN DEEP REDSHIFT-SPACE 

We analyze the structure of the two-point correlation func- 
tion in deep redshift space, using the analytic formula derived 
in the most general situation by relativistic linear perturba- 
tion theory ( Matsubara 2000 ) . In this section, the result of the 
analytic formula is briefly reviewed. Some basic notations 



of the Friedmann-Lemaitre model with dark energy extension 
are introduced in AppendixIXI 

The general two-point correlation function in redshift space 
is a function of the redshifts zi, Z2 of the two points, and 
the angle 9 between them. The analytic result of iMatsubaral 
(.2000.i) has the form, 

azuZ2,0)=D(zi)D(z2) [b(zi)b(z2)e°\zuZ2,0) 

+ f(Zl )b(Z2)e\zi , Z2 , 0) + f(Z2)b(zO^'-'Kz2 ,Zud) 

+ fiZl)fiZ2)^'-^\zuZ2,e)], (1) 

where b(z) is the redshift-dependent linear bias factor and D(z) 
is the linear growth factor normalized as D {z = ) = 1, and 
calculated by the differential equations iA5\ and ( IA6> . The 
functions are defined below. In this paper, we assume the 
concerning scales of the two-point correlation function are 
reasonably smaller than the curvature scale of the universe: 
xi2 ^ l/Tl"'/^, where jcn is the comoving distance between 
the two points. This exactly holds for a flat universe, and is 
fulfilled in practice because the curvature of the universe is 
observationally at least 20 times larger than the Hubble scale, 
c/Ho, since CIko = ^mo + ^Ao - 1 = 0.02 ± 0.02 ( Spergel et aQ 
Eooil) . This approximation does not necessarily require that 
the comoving distances to the two points from the observer, 
x(z\) and x(z2), are smaller than the curvature scale. In fact, 
the curvature effects are included in those distances in the fol- 
lowing compact version of the formula. 

Adopting the above approximation, the components ^^''^ in 
the analytic formula (Q) are given by 



e'\zuZ2,0) = C(^i2), 

, Z2 , 0) = l^i^n) + A 12 COS 7i2e;"(^i2) 

+ (^COs2 7i2-i^C2"(^12), 

(^^KzuZ2,0)= j- (l+2coH^e) d\xn) 



(2) 



(3) 
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Ai2A2iCOs6'C^'*(.Xl2) 



3 

+ [A 12 (cos 712 -2 cos 721 cos 9) 



+A21 (cos 721-2 cos 712 cos 9)] f '/'(x^) 
^ + ^ cos^ 9 - (cos^ 712 + cos^ 721 ) 



+4 cos 712 cos 721 cos9 



C^''(X12) 



+A12A2I ( COS 712 COS 721 + ^ COS 6* ) (,f\xi2) 



+ [Ai2 (5cOS7i2COS^721-COS7i2 + 2cOS72lCOs6') 

+A21 (5 COS 721 cos^ 712 - COS 721 + 2cos 712 COS 9) ] 

X^f(xi2) 



1 

+ - 
7 



1 2 

-+ - COS 

5 5 



(cos^7i2 + cos^72i) 
+4 COS 7i2 COS 721 cos9 



+7cos^7i2Cos^72i 



^f(xi2). 



(4) 
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In the above equations, the functions ^J"' are defined by 



where ji is the spherical Bessel function, and xi2 denotes the 
comoving distance between two points specified by variables 
{zi,Z2,d)- An explicit representation of xi2 is given by 

= [sHzi) + S\z2) - 2C(zi)C(z2)S(zi)S(z2)cos e 

-KSHzi)S\z2){l+cos^9)y'\ (6) 

where the first approximation is consisitent with xi2 ^ 
\K\-^'^, and 

S(z) = Sk[x(z)] 

' (-/:)-i/2sinh [(-K)^^Mz)] , (K < 0), 
[x(z), (K = 0), (7) 

^K-^/^sm[K^/h(z)], (K>0), 

,e (cosh[(-K)'/h(z)],(K<0), 
C(z)=-^(z)= <^ 1, iK = 0), (8) 

^■^ [ cos , iK>0). 

The angle 712 is the angle between the separation X[2 and the 
line of sight at the first point of zi , and 721 is the similar angle 
at the second point of Z2- Trigonometric geometry shows that 
they are analytically represented by 



cos 712 = 



S(ZI)C(Z2)-C(ZI)S(Z2)C0S9 
Sk(Xi2) 



(9) 



and cos 721 is given by a replacement of z\ ^ Z2 in the above 
equation. If the universe is exactly flat, the above equation 
simply reduces to cos 712 = (xi -X2cos0)/xi2. In this paper, 
we generally consider a non-flat universe as well as a flat uni- 
verse. The quantities A12 andA2i are needed when the redshift 
evolution of the clustering between the two points and that of 
the selection function is not negligible, and are given by 



Ai2 = SK(xi2)H(zi) — ln[H(zi)D(zi)f(ziMzi)], 
dzi 



(10) 



and the corresponding equation of A21 with a replacement 
Zi <-> Z2 in this equation. In the above equation, n(z) = dN(< 
z)/dz is the selection function including the volume factor, 
i.e., the differential number count of galaxies within a fixed 
angular area as a function of redshift z- Since H{z)d/dz = 
d jdx, the quantities A 12 and A21 have the order of relative dif- 
ference of the product H(z)D(z)fiz)n(z) between zi and Z2- 
Therefore, if the difference of this product is negligible be- 
tween two points at zi and Z2, the terms with quantities An 
and A21 can be omitted. Since the characteristic scales on 
which the factors H(z), D{z), and f{z) vary are of order of the 
Hubble scale, it is unlikely that the contributions of these fac- 
tors are important. In practice, the determination of the galaxy 
correlation function on Hubble scales is extremely difficult 
because of the high signal-to-noise ratio. A possible domi- 
nant contribution to A 12 and A21 is from the selection function 
n{z), and therefore, A12 = SK(x\2)H(z\)d\n[n{zi)]/dzi, and 
^21 = SK(xi2)H(z2)dln[n(z2)]/dz2 on sub-Hubble scales. If 
the selection function varies on scales of interest, these terms 
should be kept. When the selection function is approximately 
the same at corresponding two points at z\ and Z2, the terms 
with A12 and A21 can be omitted at all. This does not mean 
that these terms can be omitted for completely homogeneous. 



volume-limited sample, in which n{z) increases with redshift. 
Instead, those terms can be omitted when the spatial selection 
decreases with redshift and cancels the increase of the volume 
factor per unit redshift in a region of a fixed solid angle. In 
the following analysis of this paper, these terms are assumed 
to be small and omitted for simplicity. While the quantity X12 
is symmetric with the replacement of the two points 1 2, 
the quantities 712 and A 12 are not. 

The quantity corresponds to the isotropic component of 
the correlation function, since it depends only on xi2- The 
quantities ^"^'^ and ^"^^^ are relevant to distortions by the pe- 
culiar velocity field, since / — > recovers the correlation 
function in isotropic comoving space. The peculiar veloc- 
ity anisotropically distorts the correlation function even in co- 
moving space. The extended AP effect is included in the non- 
linear, anisotropic mapping from comoving space to redshift 
space (i.e., z-space) by xi2{z\,Z2,d) of equation (|6j. 

The above equations are complete set of equations we need 
to evaluate the theoretical prediction of the correlation func- 
tion in redshift space when X12 ^ l/Tl"'/^. A detailed deriva- 
tion of the equations and more general formula without the 
last restriction are given in Matsubara (2000). As discussed 
in this reference, all the known formula of the correlation 
function in redshift space are limiting cases of the general 
formula presented here. When the galaxy sample is suffi- 
ciently shallow, zi,Z2 <C 1, the formula above reduces to the 
result derived by Szalav. Matsubara & Land\J (Il998h . which 
includes wide-angle effects on the peculiar-velocity distor- 
tions. When the separation of the two-points is much smaller 
than the distances to them, X12 <C x{zi),x{z2), the formula 
reduces to the result by Matsubara & Suto ( 1996) with the 
distant-observer approximation in deep redshift space. When 
the above tw o limi ts are simultaneously applied, the result 
of HamiltoiJ (Il992h is recovered. By the Fourier transform 
of the correlation function in distant-observer approxima- 
tion, the power spectrum in deep redshift space derived by 
Ballinaer, Peacock & Heavens (1996) with distant-observer 
approximation is obtained. The last expression with shallow 
limit z ^ 1 corresponds to the original Kaiser's formula of 
the power spectrum in shallow redshift space with distant- 
observer approximation ( Kaiser. 198 7). 

3. COSMOLOGICAL INFORMATION IN THE CORRELATION 
FUNCTION IN DEEP REDSHIFT-SPACE 

The physical effects on the correlation function in deep 
redshift-space can be placed in three categories. The first one 
consists of the effects on the underlying mass power spectrum 
in comoving space. The physical parameters which primarily 
determine the shape of the power spectrum are Aj, n,, iluoh, 
/b, h, il^o, etc., where A, is the amplitude of the power spec- 
trum, n, is the primordial spectral index, h is the Hubble pa- 
rameter normalized by 100 km/s/Mpc, /b = Hbu/^mo is the 
fraction of the baryon density to the total mass density, ftiyo is 
the density parameter of neutrinos. 

The primordial spectral index n, determines the overall 
shape of the power spectrum. In the standard Harrison- 
Zel'dovich spectrum, Hs = 1 . The parameter iluoh determines 
the scale of the particle horizon at the equality epoch Zeq, since 
the radiation density is accurately known by the temperature 
of the cosmic microwave background, To = 2.725 ±0.001 K 
( Mather et al. 1999). As a result, the power spectrum has a 
characteristic peak at ^eq oc iluoh (/i~'Mpc)~'. The baryon 
fraction /b and Hubble parameter h, as well as flMoh are re- 
sponsible to the scales and strength of acoustic oscillations 
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and Silk damping. The dependences are not expressed by sim- 
ple scaling relations, a nd useful fitting relation s on physical 
grounds are provided bv lEisenstein & HulJlQgSh . The density 
parameter of neutrinos, ilj^o characterizes the free streaming 
scales by the existence of the hot dark matter. In standard cold 
dark matter scenarios, this parameter is negligible. In the anal- 
ysis of galaxy clustering, the spectral amplitude As is more 
conveniently specified by the parameter ug. The relation be- 
tween A 5 and CTg depends on other parameters explained above 
which determine the shape of the power spectrum. 

In Figure [2 the dependences of the power spectrum and 
correlation function in real space on parameters ^Moh, fs, and 
h are plotted. We consider a model with Hs = 1, ^mo = 0.3, 
/b = 0.15, h = 0.7, fii^o = as a fiducial case. Thick lines 
correspond to the fiducial model, and other lines show the 
effects of varying individual parameters, considering ^uoh, 
/b and h as independent parameters. The amplitude is nor- 
malized by fTg = 1 . The upper panels show the effects on the 
power spectrum and the lower panels show that on the corre- 
lation function, which are the (3D) Fourier transforms of the 
corresponding power spectrum. The correlation function in 
comoving space is given by the function £,^\x) of equation 

The oscillatory behavior appears in the power spectrum by 
acoustic waves before recombination epoch. Since this oscil- 
latory behavior is periodical in Fourier space, just one peak 
is appeared in the correlation function. The scale of the peak 
corresponds to the sound horizon at the recombination epoch. 
We call this peak as a baryon peak in the correlation function. 

Increasing Vty^oh shifts the peak of the power spectrum to 
the right, so that the powers on large scales are suppressed 
when (78 is fixed. Correspondingly, the correlation function 
on large scales are smaller for larger Vluoh. The scale of the 
zero-point of the correlation function decreases with this pa- 
rameter. The location of the baryon peak is changed by VLyi^h 
because the sound horizon is also dependent on this parame- 
ter 

The main effect of the parameter /b is on the strength of the 
Silk damping. Therefore, increasing /b enhances the power 
on large scales when erg is fixed. The location of the baryon 
peak is less dependent on /b, and the absolute amplitude of 
the baryon peak is predominantly dependent on /b. The zero- 
point of the correlation function is not much affected as long 
as/B^O. 

Although the effect of h on the power spectrum does not 
seem to be significant in the Figure, the phase of the baryon 
oscillation is shifted. This is partly because we fix Vluoh in- 
stead of the physical dark matter density iluoh^- Since the 
fraction /b is also fixed, the physical densities varies with h 
as Omo/i^ (X h, Vl^oh^ oc h in our choice of the independent 
parameters. Although the shape of the power spectrum is de- 
termined by the physical density parameters, the length scale 
should be measured in units of /!"'Mpc in redshift space, and 
one can not eliminate the explicit dependence on the Hubble 
parameter by scalings. It is only when baryons are absent that 
the shape of the power spectrum is characterized by a single 
parameter 0,uQh, without any explicit dependence on h. 

The effect of h on the correlation function is more notice- 
able since the multiple baryon wiggles in the power spectrum 
cumulatively contribute to the baryon peak in the correlation 
function. The sound horizon at the recombination epoch in- 
creases with h when iluoh and /b are fixed, and the location 
of the baryon peak is shifted to the right. 



The second category of the physical effects on the correla- 
tion function in deep redshift-space is the velocity distortions, 
or the Kaiser's effects. These effects are included in second 
and third terms of equation Q. In linear regime, the coherent 
infall toward centers of density maxima flatten the clustering 
pattern along the lines-of-sight. Depending on the geometri- 
cal angles among the observer and the two points, the corre- 
lation function in redshift space is a linear combination of the 
functions ifixn). The functions ^'^ and are plot- 
ted in Figure|2]with the fiducial model of the power spectrum 
described above. Other functions, ^{'\ ^j^', and ^3'' are 
not significant when A 12, A21 are small and are omitted in the 
Figure. All the functions ^J"* have irregularities on the scale 
of the baryon peak. Since the correlation function in redshift 
space is expressed by a linear combination of these functions, 
the resulting function has the same irregularities on the same 
scale. 

The third category of the physical effects on the correla- 
tion function in deep redshift space is the extended AP ef- 
fect. The irregularities on the scale of baryon peak play 
a role of spherical objects with known radius in comoving 
space. Originally, Alcock & Paczvnski ( 1979) proposed that 
measuring the ellipticity of some objects which are spheri- 
cal in comoving space, the value of the cosmological constant 
can b e constrained. It was pointed o ut that the power spec- 
trum (Ballinger. Peacock & Heavens 1996) and the correla- 
tion function (Matsubara & Suto 1996) can serve as such ob- 
jects. The baryon wiggles in the power spectrum can be used 
for th is purpose ( Blake & Glazebrook.2003: Seo & Eis ensteiiil 
l2003h . While the baryonic feature is split into many wiggles 
in the power spectrum, there is just one peak in the correlation 
function. 

In Figure |3] contour plots of the resulting correlation func- 
tion in redshift space with z\ =0.1, 0.3, 1.0, and 3.0 are 
shown. The fiducial model parameters with fiKO = 0, w = -1, 
i^MO = 0.7, «.5 = 1, 57^0 = 0, f^Mo/! = 0.21, /b =0.15 are adopted. 
The bias parameters are assumed to be b{z) =1, 2, 3, and 
3 for z = 0.1, 0.3, 1.0, and 3.0, respectively. The normal- 
ization of the power spectrum is chosen so that the galaxy 
normalizations at corresponding redshifts are unity; a^giz) = 
D(z)b{z)a'i = 1. First point labeled by zi is sitting at the cen- 
ter of each plot. The contour lines represent the value of the 
correlation function in redshift space, S,{z\,Z2,0), and the co- 
ordinates of the plots are (z2sin(?,Z2Cos0). In other words, the 
plots show the contours of the correlation function in apparent 
z-space. A prominent feature in these plots is the existence of 
ridges. These ridges are perfectly circular in comoving space, 
i.e., they are located at the lines of constant X12. The scales of 
the baryon peak in comoving correlation function, which we 
have seen in the Figure [2 is indicated by thick dotted lines 
in the plots. Obviously, the ridges correspond to the baryon 
peaks in comoving correlation function. Therefore, we will 
call those ridges as the baryon ridges. While the peculiar ve- 
locities alter the amplitudes of the correlation function along 
the ridges depending on the direction relative to the lines of 
sight, the shape of the ridges are not distorted by them. There- 
fore, the baryon ridges are ideal statistically spherical objects 
which are useful for the extended AP test. The Fourier coun- 
terpart of the baryon ridge is the a coustic rings in the power 
spectrum in redshift space (iHu & H aiman 2003). While the 
rings in the power spectrum are spread over many scales, the 
baryon ridges in the correlation function have just a single 
scale. 
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Fig. 1. — The power spectram (upper panels) and the correlation function (lower panels) in comoving space. Thick lines correspond to model parameters 
Quo = 0.2, /b = 0.15, and h = 0.7. The amplitude is normalized by erg = 1. Thin lines show the variations of each parameters as suggested in the plots. When 
SImo'' increases, the amplitude on large scales decreases. When /b increases, the amplitude on large scales increases. When h increases, the location of the 
baryon peak is shifted to large scale. 
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Fig. 2. — Examples of the function 5," (x) that is needed in the analytic 
formula of the correlation function in redshift space. 



Besides the linear power spectrum at present, the correla- 
tion function in deep redshift space depends on other cos- 
mological parameters such as 51ko, ^mo, and w. The overall 



amplitude of the correlation function depends on the linear 
growth rate D(z). The peculiar velocity effects depend on the 
logarithmic growth rate f{z). The AP effects depends on the 
time-dependent Hubble parameter H{z) and comoving angu- 
lar diameter distance, Da{z) = S{z)- The function H(z) spec- 
ifies the scales along the lines of sight in redshift space, and 
the function D\(z) determines the scales perpendicular to the 
lines of sight in redshift space. Figure |4] shows the depen- 
dences of the relevant functions on cosmological parameters. 
In the derivatives by iiKo, the mass density parameter iluo is 
fixed so that the derivatives with respect to 57k() is equivalent 
to the derivatives with respect to the dark energy density pa- 
rameter, r^Qo. The redshift evolution of dark energy equation 
of state is parameterized by w = wo + wiz/(l+z)- The parame- 
ter sensitivity of the physical quantities are generically higher 
in deep redshift samples. The dark energy is less important 
when the redshift is too high. Only the function f(z) is sen- 
sitive to i7M0 at low redshifts. Overall, parameter sensitivity 
is high at around z ^ 1.0. Therefore, intense redshift surveys 
around redshift ^ 1 provide important constraints on cosmo- 
logical models. 

4. CONSTRAINING COSMOLOGICAL PARAMETERS WITH DEEP 
REDSHIFT SAMPLES 

4. 1. Evaluating the Fisher Information Matrix 

Accurate estimations of how the correlation function in 
deep redshift space can constrain the cosmological parame- 
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Fig. 3. — Contour plots of the coiTelation function l;(zi,Z2,S) in deep redshift space. Tlie first point with redshift zi is sitting at the center of the plot. The 
contours represent the correlation value as a function of za and S and the coordinates are defined by {z2 sin 0, Z2 cos 8). In other words, the plots show the contours 
of the correlation function in apparent z-space. The fiducial model defined in the text is assumed. 



ters are crucial in designing redshift surveys. For this purpose, 
we assume idealized surveys which are accessible by present- 
day technologies, and give estimates of the error bounds for 
cosmological parameters, calculating Fisher information ma- 
trices. In this section, we consider JIko, ^^mo, ^Moh, /b, 
w, and b as an independent set of cosmological parame- 
ters. The bias is the most uncertain factor in galaxy redshift 
surveys. On linear scales, the complex bias uncertainty is 
renormalized to a scale-independent linear bias parameter, b 
(|§cherrer & Weinberg 1998; Matsubara 1999). The bias pa- 
rameter is not a universal parameter and varies from sample 
to sample, depending on selection criteria of individual ob- 
servations and also on redshifts. In the following, we always 
consider the situation that the bias parameter in a given sample 
is unknown in each sample and should be determined simul- 
taneously with other parameters. In other words, the bias is 
always marginalized over in the Fisher analysis below. There- 
fore, the results in the following are free from the bias uncer- 
tainty unless the variation of the bias in each sample is too 



strong. 

The forecasts of expected errors in parameter determina- 
tions in a given sample are most easily obtained b y the anal- 
ysis of the Fisher info rmation matrix (e.g., Kendall & Stuarll 
fl969H TTi errierillT992h . This technique enables us to obtain 
expected errors in parameter determination by a given sam- 
ple without performing any Monte-Carlo simulation. The 
Fisher information matrix is the expected curvature matrix 
around a maximum point of the logarithmic likelihood func- 
tion £ in parameter space as defined by the equation iBli . 
When the distribution of the data is given approximately 
by a multivariate Gaussian one wit h a correlation matrix C, 
the Fisher matrix is given by (e.g ., FVogelev & Szalavlfi99d: 
iTe gmark. Tavlor & Heavensll997h 



= ^Tr ( C 



dC 



-C" 



dC 



(11) 



where {6 a} is the set of model parameters. Predictions of the 
error bounds when some parameters are marginalized over are 
also useful. For this purpose, we can define the marginalized 



7 



0.5 



-0.5 



0.4 



a; 0.2 



-0.2 



-'"1 1 1 - 

dlnH(z)/dQj„ 

- dln£l,(z)/dn^„ 

dlnZ)(z)/dn,„ 

dln/(z)/dn,„ 

X 


dlnH(z)/dlnn„„ - 

din D,(z )/dln n„„ - 

dlnD(z)/dln n„„ 

\ dln/(z)/dln n„„ \ 


■ ■.. \ ..■'■\ ■ 

: ■■■■ \^ 




'w\ \ — 1 1 1 1 III \ — 1 1 hill 

dlnHCzj/difj 

dlni),(z)/dw„ 

dlnC(z)/dvir„ 

dln/lzj/dWj 


'W\ 1 — 1 1 1 III 1 — 1 1 1 1 II 

dlnHizj/dw, 

dlnD,(z)/dWi 

dlnD(z)/dw, 

dlnf(z)/dw, 


^ ~^ / — " " " 
^ x ^ ^ C- ^ ^ 

M.I 1 . TiTi'l 1 - 


M,l 1 1 - 



0.1 



10 0.1 



10 



Fig. 4. — Sensitivities of the cosmological functions on parameters. To present the fractional sensitivities, derivatives of the logarithm of the Hubble parameter 
H{z) (dotted lines), the comoving angular diameter distance Da(z) (short dashed lines), the linear growth rate D(z) (long dashed lines), and logarithmic growth 
rate /(z) (dot-dashed lines) are shown as functions of the redshift. Upper left panel shows the sensitivities on the curvature parameter f2Ko> fixing f!MO- Upper 
right panel on the density parameter fiMO. fixing Oko- Lower panels on the parameters of equation of state of the dark energy, wq and w\ . 



Fisher matrix by a method presented in AppendixlEl and ob- 
tain the corresponding error forecasts. The bias parameters 
are always marginaHzed over by that method. 

Thus, the calculation of the Fisher matrix is straightforward 
once the correlation matrix C of the data and its derivatives 
are obtained. However, the calculation of the correlation ma- 
trix depends on what kind of observed data is analyzed. A 
straightforward definition of the data in an analysis of the two- 
point correlation function is the binned values of the correla- 
tion function. The two-point correlation function in redshift 
space is a function of a set of the variables {z\^zi.,&)- With 
a suitable binning of this three-dimensional space, we obtain 
the data vector ^ (assumed to be a column vector), which con- 
sists of the values of the correlation function averaged over 
within each bin. In this case, the correlation matrix is given 
by (7 = (^^^). To evaluate this correlation matrix, appropriate 
estimates of covariances between components of £ in a given 
observation with certain geometry and selection function are 
necessary. This is marginally possible, but requires intensive 
multidimensional integrations to obtain the proper correlation 
matrix. 

Fortunately, there is a simpler, alternative method of 
evaluating the Fisher matrix of a given sample, us- 
ing the information of the two- point correlation function 
jMatsubara & Szalavll200l 120031) . In this method, the data 
vector is taken to be pixelized galaxy counts A^,- in a survey 
sample. Thus, the correlation matrix is simply given by the 
smoothed correlation function convolved by pixels: 

C,v = ((M-^0(^;-^>)) 



d\i f dhi azi,z2,en)+mj, (12) 

where v, is the volume of the pixel / and .^V, is the ex- 
pected number of galaxies in that pixel. The integrated re- 
gions are limited within pixels in redshift space. Th e sec- 
ond term is the contribution from the shot noise ( Peeblej 
1980). At first sight, the six-dimensional integration in equa- 
tion ( I12t se ems computationally co stly to perform. Never- 
theless, Matsubara. Szalav. & Popel 42004) showed that there 
is no need to perform these direct integrations when one can 
set the spherical shape of the pixels in comoving space. In 
which case, performing the smoothing integrations of equa- 
tion ( I12t is equivalent to just replacing the power spectrum 
P(k) by P(k)W^(kR) in the equation Q, where W(x) = (sinx- 
xcosx)/3x^ is the smoothing kernel in Fourier space and R is 
the smoothing radius in comoving space. On one hand, one 
need to know the cosmological parameters i^MO, ^qo, and w 
in advance to set pixels in redshift space which are exactly 
spherical in comoving space. On the other hand, the slight el- 
lipticity of the pixels does not change the value of the correla- 
tion matrix which is calculated by assuming spherical pixels. 
Therefore, with an approximate set of cosmological parame- 
ters, one need not to perform the integrations of equation il2\ 
and directly obtain an accurate correlation matrix in just the 
same way as computing the correlation function itself in red- 
shift space. 

In the following, we use the latter method to calculate the 
Fisher matrix. The correlation matrix of equation il2\ is cal- 
culated from the expected number density in each cells and 
models of the correlation function. The Fisher matrix with 
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Fig. 5.— Illustration of an optimal survey design with multi-layers of pm g _ Another survey design with continuous layers, 

regions. 



any marginalization is obtained only from the correlation ma- 
trix. 

4.2. Baseline Redshift Surveys 

One of the advantages of the deep redshift surveys is that 
the information from the time sequence of the galaxy cluster- 
ing is available. Since the solid angle on the sky is limited 
to 47r at most, the comoving volume in which we can observe 
galaxies is larger in high-redshift universe than in low-redshift 
universe. If the galaxies in a spectroscopic survey are selected 
in a fixed solid angle by certain criteria, as commonly done 
in redshift surveys like the 2dF and the SDSS surveys, the 
number density of observed galaxies becomes sparse at high 
redshifts. This is a critical drawback in a correlation anal- 
ysis since the shot noise dominates the clustering signal on 
scales of interest. For example, SDSS quasar sample is not 
optimal to determine the cosmological parameters from the 
clustering, because of their sparseness ("Matsubara & Szalav 
12002). To maximally take advantage of the linear formula 
of the correlation function in deep redshift space, clustering 
properties on scales of 10-200 /i~'Mpc provides a useful in- 
formation on cosmology. Too small scales suffer the nonlin- 
ear effects which has less cosmological information and have 
difficulties in analytical treatments. Therefore, observations 
with a number density of over about 1 /(10/;~'Mpc)'' is ideal. 

To retain comparable number densities at different red- 
shifts, the survey region on the sky should be narrower at 
high redshifts than at low redshifts, since the total number of 
observed galaxies is usually limited by the observation time, 
or budgets. Therefore, an optimal design of the redshift sur- 
vey will be something similar to that illustrated in Figure |5l 
The optimal survey illustrated here consists of several layers 
of sub-surveys. Each sub-survey is optimized to catalog par- 
ticular redshift range and has selection criteria that is most 
effective to pick the galaxies up in corresponding range of 
redshift. While simple color selections can be used for this 
purpose, the photometric redshift data might be the best to ful- 
fill the required selections. In low-redshift layers, the surface 
number density of galaxies for the spectroscopy can be low, 
while the sky coverage of the survey regions should be wide. 
In high-redshift layers, the surface number density should be 
high, while the sky coverage can be narrow. The requirements 
of the telescope for each layers are different. Thus this opti- 
mal survey would effectively carried out by several telescopes 
with a wide spectrum of capabilities such as resolutions and 
wideness of the field of view. 

One can also imagine a limit of infinite number of layers. 
The survey geometry in this case is like a big column in red- 
shift space, as illustrated in Figure |6l In this case, the selec- 
tion criteria should continuously vary with the position on the 
sky, therefore delicate selections are required to sample galax- 
ies as homogeneously as possible. While technically more 



challenging, this survey strategy is ideal for the cosmologi- 
cal analysis of the clustering at various redshifts in the deep 
universe. 

In this section, however, we consider more realistic surveys 
that can be achieved by present-day technologies. We take 
the series of baseline redshift surveys explained below. The 
surveys are categorized to four types, which are assumed to 
be volume-limited samples. The choices of samples at high- 
redshifts are similar to that considered by .Seo & Eisensteii] 
(2003). 

The first sample is a low-redshift sample around z ~ 0.15. 
This sample is motivated by a volume-limited su bsample of 
galaxies in the SDSS survey (Teemark et al. 2004) with ab- 
solute r-magnitude M,- ~ -22. Our baseline sample has a 
uniform number density n = 1 x 10~^/(/!~'Mpc)^ in the red- 
shift range z = 0.05-0.25. The sky coverage is assumed to be 
f^sky = 8.05 str. The bias factor is assumed to be = 1.3. We 
call this sample as Sample A below. 

The second sample is the one aroundz ~ 0.3. This sample is 
motivated by the volume-limited catalog of the Luminous Re d 
Galaxies (LRGs) in the SDSS survey iEisenstein et alJ200ll) . 
Our sample is assumed to have n = 1 x \Qr^ / {hT^yi^cf , z = 
0.2-0.4, Osky = 2.57 str, and b = 2.0. We call this sample as 
Sample B below. 

The next samples are the series of surveys around z ^ I- 
The choice of the target galaxies is not t rivial. We follow 
the ba seUne surveys that are considered bv lSeo & Eisenstei^ 
J2003h . The target galaxies are assumed to be either gi- 
ant ellipticals or luminous star-forming galaxies. This cat- 
egory of samples consists of surveys of four redshift bins, 
z = 0.5-0.7, z = 0.7-0.9, z = 0.9-1.1, and z = 1.1-1.3. The 
sky coverage is assumed to be ilsky = 0.914, 0.647, 0.517, and 
0.451, respectively. The bias factors are 1.25, 1.40, 1.55, and 
1 .70, respectively. Approximately we have crgg ~ 1 in each 
samples with the above choices of bias factors. We assume 
n = 5 X lQ~'^/{h~^Mpc)^ commonly for each samples. We call 
the series of samples as Sample CI to Sample C4 below. 

The last sample is the survey of Lyman break galaxies 
around z ^ 3. The selection techniques for this kind of galax- 
ies are reported bv lSteidel et alJ ( 1 19961) . We consider the sam- 
ple with a redshift range of z = 2.65-3.35. We assume b='i3 
and n=\ x 10~^/(/?~'Mpc)^ which are consistent w ith the ob- 
servation ('SteideLitaIIl998*'Ade^^ We call 
this sample as Sample D. The parameters of the baseline sur- 
veys are summarized in Table [2 

The Fisher matrix is calculated by assuming a cone geom- 
etry with an open angle which gives 1/10 of the assumed sky 
coverage in each sample. The spherical cells with a comov- 
ing radius 15/;~'Mpc are placed in closed-packed structure in 
a survey region. The total number of cells is exactly 4900 
in each 1/10 subsample. Thus the comoving volumes of all 
samples are essentially the same. The resulting Fisher ma- 
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TABLE 1 

Parameters of the baseline surveys. 



Samples 


z. 


area [str.] 


h 


n [/(lOr'MpcV'] 


Candidates 


A 


0.05-0.25 


8.05 


1.30 


1.0 


M,- ~ -22 galaxies 


B 


0.2-0.4 


2.57 


2.00 


0.1 


Luminous red galaxies 


CI 


0.5-0.7 


0.914 


1.25 


0.5 




C2 


0.7-0.9 


0.647 


1.40 


0.5 


Giant ellipticals, or 


C3 


0.9-1.1 


0.517 


1.55 


0.5 


star-forming galaxies 


C4 


1.1-1.3 


0.451 


1.70 


0.5 




D 


2.65-3.35 


0.102 


3.30 


1.0 


Lyman break galaxies 



trix is multiplied by 10 to obtain the final Fisher matrix. This 
method of subdividing the sample into 1/10 is taken to re- 
duce the dimension of the matrix and to compromise with the 
CPU time and the memory requirement of handling the large 
correlation matrix. In general, the subdivision overestimates 
the error bounds because the information from cross correla- 
tions between subsamples is neglected. However, the larger 
the sample is, the less such cross correlations are. We con- 
firmed that the error bounds of the full sample is accurately 
calculated by just adding the Fisher matrices of the subsam- 
ples in our case. 

4.3. Results of the Fisher Analysis 

According to the calculation of the Fisher matrix, the ex- 
pected error bounds with any marginalization are straightfor- 
wardly obtained by methods explained in AppendixIBI Since 
the bias parameter should be simultaneously determined sam- 
ple by sample, we calculate the Fisher matrix of each sample 
with bias separately marginalized over within each sample. 
Therefore, the results below are not affected by bias uncer- 
tainties, as long as the bias factor does not significantly varies 
within each samples. 

In Figure concentration ellipsoids in 2-parameter esti- 
mations are plotted, i.e., each panel shows the expected er- 
ror bounds when corresponding two parameters are simulta- 
neously determined, fixing other parameters except the bias 
parameter. The bias parameters are marginalized over. The 
Fisher matrix of the combined sample are also calculated by 
adding the Fisher matrices of Sample A-D with bias sepa- 
rately marginalized over sample by sample. 

Most of the pairs of parameters do not exhibit degener- 
acy between them. This means that those pairs of parame- 
ters contribute to the correlation function fairly differently as 
explained in §3. However, one can notice relatively clear de- 
generacies between iiKO and w. Since we consider r^MO and 
r^Ko as independent parameters, changing the parameter JIko 
is equivalent to changing the dark energy density figo with 
f^MO fixed. The parameters of the dark energy are mainly con- 
strained by AP effect besides the linear growth rate. In the 
low-redshift samples, the parameters J7qo and w only depend 
on the growth rate and thus degenerate with each other Since 
the bias parameter and the growth rate is distinguishable by 
the velocity distortion, one can still extract information on the 
dark energy from low-redshift sample with the cost of the de- 
generacy. On the other hand, high-redshift samples can actu- 
ally constrain these two dark energy parameter independently. 

The direction of the major axis of the concentration ellipse 
in the VLyj^-w plane rotates anti-clockwise with the average 
redshift of the samples. Accordingly, the concentration el- 



lipse of the combined sample is quite small. For example, the 
error bound of VIqq is ~ ±2% and that of w is ^ ±5%. The 
presented error bounds are given when the other parameters 
are fixed except the bias. It is also interesting to see the ex- 
pected constraints when all parameters are determined by only 
using the deep galaxy surveys. In Figure|8] the concentration 
ellipses are plotted on VlYia-w plane when the other parame- 
ters iluo, ^Moh, /b, cg and b's are all marginalized over The 
complementarity of the surveys at various redshifts is obvious 
in this case. The concentration ellipses of high-redshift sam- 
ples are also elongated by uncertainties of other parameters. 
The direction of the major axis of the ellipses rotates anti- 
clockwise with redshifts. Consequently, the degeneracy be- 
tween i7K0 and w are broken in the combined sample and the 
error bound is much smaller than in the individual samples. 
The error bound of ftqa is ~ ±4% and that of w is ~ ±10% 
even when all the other parameters are completely unknown. 

5. CONCLUSIONS AND DISCUSSION 

In this paper, properties of the two-point correlation func- 
tion in deep redshift space is theoretically investigated. A 
completely gene ral expression of th e linear correlation func- 
tion derived by 'Ma tsubaral iZOOdh is reduced to a simpler 
form, ea.(lll- (llO> with an approximation which is valid for 
realistic surveys. Baryon wiggles in the power spectrum cor- 
respond to a single peak in the correlation function. This peak 
differently depends on the parameters of the underlying power 
spectrum, ^M{)h, /b, and h. In the two-dimensional contour 
plot of the correlation function in redshift space, the corre- 
sponding feature is the baryonic ridge. The shape of the ridge 
is perfectly circular in comoving space, in spite of the effects 
of the peculiar velocity. This is preferable for the cosmolog- 
ical test by the AP effect. Using the analysis of the Fisher 
information matrix, the cosmological test by the correlation 
function in deep redshift space is shown to constrain the prop- 
erties of the dark energy component as well as other cosmo- 
logical parameters, even if the bias is uncertain. In particular, 
there is a clear complementarity of the samples of various red- 
shift ranges in probing the dark energy component. 

Therefore, a survey of a wide redshift range is particularly 
useful to probe the dark energy. Keeping a sufficient number 
density of galaxies in the given observation time, examples of 
the optimally designed survey geometry are illustrated in Fig- 
ure |5] and Figure |6l Selecting galaxies for the spectroscopic 
follow-ups as densely as possible is one of the major techni- 
cal challenges in this kind of survey strategy. However, recent 
advances of the technology enable us to carry it out. With 
the advent of 8-lOm telescopes, galaxies of 1.4 can be 
straightforwardly selected by color selections as demonstrated 
by, e.g., the SDSS survey jYork et a l."2000':'Eisenstein & Hd 
1998) and the DEEP2 Galaxy Redshift Survey (Davis et aj 
2003; Coil et al. 2003). In z ^ 2.5, Lyman-break galax ies are 
photometrically selected with sufficient density ( Steidel et alJ 
ll996lfT99llAdel berger et al. 1998). The range 1.4 < z < 2.5 
is known as the redshift desert, since spectroscopic iden- 
tifications o f galaxies were historically difficult. Recently, 
ISteidel etal . (2004) showed that large numbers of galaxies in 
this range can actually be selected with standard broadband 
color selection technique using the Keck I Telescope with the 
LRIS-B spectrograph. Thus galaxies at z ^ 3 can be opti- 
cally selected by ground-based 8- 10m telescopes with suffi- 
cient number densities ( Adelberger et al. 2004). Therefore, 
the proposed surveys illustrated by Figure |5] and Figure |6l in 
the redshift range of z ^ 3 with approximately homogeneous 
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Fig. 7. — Two-parameter joint error bounds expected by the baseline surveys. Each ellipse represents the expected constraint on the parameter space with 
1(7 significance. In each panel, corresponding two parameters are varied with bias marginaHzed over Other parameters are fixed in each plot. Error bounds for 
Sample A (dotted lines), Sample B (short-dashed Unes), Sample C1-C4 (long-dashed lines). Sample D (dot-dashed lines) are plotted together with a combined 
sample (solid lines). 



selections are accessible by present-day technologies. 

Cosmology with the large-scale structure has mainly been 
studied by the snapshot of the nearly present universe, since 
the sufficient number of galaxies have been available only at 
shallow redshifts, z ^ 0.3. CMB observations, on the other 
hand, map the universe of the z— 11 00 universe. The obser- 
vations of the galaxy clustering in the deep-redshift universe 
of 0.3 ^ z ^ 3, which are within reach of the present-day tech- 
nology, will provide a unique information on the evolving uni- 
verse with the analysis of the correlation function presented in 
this paper. Observations of the evolving universe is momen- 
tous not only to probe the nature of the dark energy compo- 
nent of the universe, but also to establish or even falsify the 
standard picture of the cosmology. Only data of the snap- 



shots of the universe are not enough to confirm the consis- 
tency of the cosmological model. Whether or not the structure 
of the universe with respect to both space and time are fuUy 
explained by minimal set of cosmological assumptions and 
parameters will be a major concern in cosmology. In the era 
of the precision cosmology, any need for additional assump- 
tions will suggest terra incognita in the physical universe. 



1 thank Alexander Szalay and Adrian Pope for helpful dis- 
cussion. 1 acknowledge support from the Ministry of Educa- 
tion, Culture, Sports, Science, and Technology, Grant-in- Aid 
for Encouragement of Young Scientists, 15740151, 2003. 



11 



-0.5 



-1.5 



1 1 ''!'''!/' 
/' 

/ 

/ 

/ 


V / / // y / 
A/ / // / /.-^- 


/' 

/' 

/■ / 

- 


/y/lJ'/ y ' 
///,//7/.-/ - 

H-i// 

" ' //// f^" 




\ // y 

fy .y 


- A-fO / 


r : 

/' 


ijyy/(\ , 





-0.4 



-0.2 



0.2 



0.4 



Fig. 8. — Joint error bounds of the dark energy parameters expected by the 
basehne surveys. Parameters Cmo. f^Mo''. /b, fg and b are all marginalized 
over. The correspondences between lines and samples are the same as in 
FigurelTl 

APPENDIX 
A. BASIC NOTATIONS 

In this Appendix, basic notations of the Friedmann-Lemaitre universe with dark energy extension used in this paper are intro- 
duced. 

The background metric is given by the Robertson- Walker metric, 

ds^=-dt^ + a\t) [dx^ + SK\x){de^ + sm^e(t?)] , (Al) 
where we employ a unit system with c = 1, and adopt a notation, 

{-K)-^l\m\\ [i-Kyi^x] , (K < 0), 

X, iK = Q), (A2) 

sin [/T'/^x] , (^:>0), 

where K is the spatial curvature of the universe. The comoving distance x(z) at redshift z is given by 

dz! 



Sk(x) ; 



Hiz'y 



(A3) 



where H(z) is the redshift-dependent Hubble parameter When we allow the dark energy component having a non-trivial equation 
of state, p(z) = w(z)piz), the Hubble parameter is given by 



H(z) = Ho 



(i+zfnMo-(i+zfnKQ+sxp{3 



1+W 



dz]n 



QO 



1/2 



(A4) 



where f^MO is the density parameter of matter, JIqo is the density parameter of dark energy, and ilKO = ^^mo + ^^qo - 1 is the 
curvature parameter. The linear growth factor D{z) and its lo garithmic derivative, f(z) = dlnD/dlna = -(1 +z)dlnD/dz are the 
solution of the following simultaneous differential equations iMatsubara & Szalavl2003h 

dlnD 

= /, (A5) 



d\na 
df . 
din a 



l + 3w 



— — r!Q /+-r!M, 



(A6) 



where the scale factor a is related to the redshift by a = (1 +z) and J^m, ^^q are the time-dependent density parameters of matter 
and dark energy, respectively: 



^m(z)= 77^(1 +z)'f^MO 



^q{z) = 



Ho' 
HHz) 



exp I 3 



1 + w 

( 

1+z 



dz riqo 



(A7) 
(A8) 



The normalization D(z = 0) = 1 is adopted throughout this paper. 
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B. MARGINALIZED FISHER MATRIX 

In this Appendix, useful equations to obtain estimates of the error covariance from the Fisher matrix when some parameters 
are marginalized over The Fisher matrix F is an expectation value of the curvature matrix of the logarithmic likelihood function 
in parameter space, 

where = {0\,02, ■ ■ .)^ is a set of model parameters, and C is the likelihood function. In the Fisher matrix analysis, the error 
covariance matrix (A0A0^) of a set of model parameters asymptotically corresponds to the inverse of the Fisher matrix 

F"' ~ (A6»A6»^), (B2) 

when all the parameters are simultaneously determined. 

The error covariance matrix of partial set of parameters, marginalizing over other parameters is given by a sub-matrix of the 
inverse of the Fisher matrix. Therefore, a Fisher matrix F of these partial set of parameters corresponds to the inverse of the 
sub-matrix of equation ( IB2> . We consider the situation that the first n parameters are estimated and other m parameters are 
marginaUzed over The full Fisher matrix has a form. 



where A is an n x n symmetric matrix, B is an n x m matrix, and C is an m x m symmetric matrix. A representation of the inverse 
of the equation JB4l i is given by 

^-1 _ I {a-bc-^bY -{a-bc-^bYbc-^ \ 

\^-(C-BTa-'B)"'bTa-i (C-BTa-'b)-' )' 

as can be explicitly confirmed that F~^F = I. Therefore, the marginalized nx n Fisher matrix F of the reduced parameters 
corresponds to 

F =A-BC~^B^ . (B5) 

The first term is the Fisher matrix of the reduced parameter space, without marginalization, and the second term corresponds to 
the effects of the marginalization. 

Given a Fisher matrix F, with or without marginalization, the estimations of the error bounds are obtained by contours of 
concentration ellipsoid /S.9^F/S.9 = const., where const, depends on significance levels and dimensionality of the parameter 
space. Therefore, obtaining the marginalized Fisher matrix by the equation iB5i is sufficient to obtain the error estimation by 
drawing concentration ellipsoids, and calculation of an inverse of the full Fisher matrix is not necessary. 

Below, we derive the correspondence between contour levels of AO^FAO and the significance level in the joint estimation of 
multiple parameters. Assuming a Gaussian profile around a peak of the likelihood function, and the correspondence between 
covariance matrix and the Fisher matrix of equation iB2\ . an estimate of j/cr error bounds can be obtained by solving an equation 
foTX(iy), 

(27r)"detF-' / d"9exp{--Ae'^FAe)=^= g-''/^ (B6) 

where _ 

Ae^FAe=X^(y) (B7) 

defines a concentration ellipsoid of vg significance level. Adopting a linear transformation of variables, x = L7A9, where L is a 
nxn matrix by a Cholesky decomposition of the Fisher matrix, F = LL7 , the equation ( IB6> reduces to a simple equation, 

l-X 



t/xx"-ie--''/2 = 2"/2-ir(n/2)erf(i//\/2) , (B8) 



where erf(x) = 2tt '/^ J^dte '' is the error function. Numerical solutions of this equation in several cases are given in Table IB2 
In this paper, la concentration ellipses in marginalized 2-parameter space are presented, i.e., n = 2,v=\ and X = \ .52. 

REFERENCES 

Adelberger, K. L., Steidel, C. C, Giavalisco, M., Dickinson, M., Pettini, M., Coil, A. L., et al. 2003, ApJ, in press (ArXiv pto)-ph/0305586l 

&Kellogg, M. 1998, ApJ, 505, 18 Colless, M. et al. 2001, MNRAS, 328, 1039 

Adelberger, K. L., Steidel, C. C, Shapley, A. E., Hunt, M. R, Erb, D. K., Davis, M., et al. 2003, Proc. SPIE, 4834, 161 (ArXiv |astto^ph/0209419> 

Reddy, N. A., & Pettini, M. 2004, ApJ, 607, 226 Eisenstein, D. J. & Hu, W. 1998, ApJ, 496, 605 

Alcock, C. & Paczynski, B. 1979, Nature. 281, 358 Eisenstein, D. J. et al. 2001, AJ, 122, 2267 

Ballinger, W. E. & Peacock, J. A. & Heavens, A. E 1996, MNRAS, 282. 877 Hamilton, A. J. S. 1992, ApJ, 385, L5 

Bennett, C. L. et al. 2003, ApJS, 148, 1 Hu, W. & Haiman, Z. 2003, Phys. Rev. D, 68, 063004 

Blake, C. & Glazebrook, K. 2003, ApJ, 594, 665 Kaiser, N. 1987, MNRAS, 227, 1 



13 



TABLE B2 

Contour levels X{i/) of the concentration ellipsoids in joint estimations of parameters. 





n= 1 


n = 2 


n = 3 


I/ = 1 


1 


1.52 


2.17 


i/ = 2 


2 


2.49 


3.12 


u = 3 


3 


3.44 


4.03 



Kendall, M. G. & Stuart, A. 1969, The Advanced Theory of Statistics, Vol. 2 

(London: Griffin) 
Linder, E. V. 2003, Phys. Rev. D, 68, 083504 

Mather, J. C, Fixsen, D. J., Shafer, R. A., Mosier, C, & Wilkinson, D. T. 

1999, ApJ, 512, 511 
Matsubara, T. 1999, ApJ, 525, 543 
Matsubara, T. 2000, ApJ, 535, 1 

Matsubara, T, Szalay, A. S., & Pope, A. C. 2004, ApJ, in press (ArXiv: 

astro-ph/0401251 i 
Matsubara, T. & Suto, Y. 1996, ApJ, 470, LI 
Matsubara, T. & Szalay, A. S. 2002, ApJ, 574, 1 

Matsubara, T. & Szalay, A. S. 2003, Physical Review Letters, 90, 021302 
Peebles, P. J. E. 1980, The Large-Scale Structure of the Universe (Princeton; 

Princeton University Press) 
Perlmutter. S. et al. 1999, ApJ, 517, 565 
Riess, A. G. et al. 1998, AJ, 116, 1009 
Saunders, W. et al. 2000, MNRAS, 317, 55 
Scherrer, R. J. & Weinberg, D. H. 1998, ApJ, 504, 607 
Seo, H. & Eisenstein, D. J. 2003, ApJ, 598, 720 



Shectman, S. A., Landy, S. D., Oemler, A., Tucker, D. L., Lin, H., Kirshner, 

R. P, & Schechter, P L. 1996, ApJ, 470, 172 
Spergel, D. N. et al. 2003, ApJS, 148, 175 

Steidel, C. C, Giavalisco, M., Pettini, M., Dickinson, M., & Adelberger, K. L. 

1996, ApJ, 462, L17 
Steidel, C. C, Adelberger, K. L., Dickinson, M., Giavalisco, M., Pettini. M., 

& Kellogg, M. 1998, ApJ, 492, 428 
Steidel, C. C, Shapley, A. E., Pettini, M., Adelberger, K. L., Erb, 

D. K., Reddy, N. A., & Hunt, M. P 2004, ApJ, 604, 534 (ArXiv; 

astro-ph/0401439i 
Szalay, A. S., Matsubara, T, & Landy, S. D. 1998, ApJ, 498, LI 
Tegmark, M., Taylor, A. N., & Heavens, A. F. 1997, ApJ, 480, 22 
Tegmark, M., et al. 2004, ApJ, 606, 702 

Therrien, C. W. 1992, Discrete Random Signals and Statistical Signal 

Processing, (New Jersey: Prentice-Hall). 
Vogeley, M. S. & Szalay, A. S. 1996, ApJ, 465, 34 
York, D. G. et al. 2000, AJ, 120, 1579 



